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It is well known that, when an external general relativistic (electric-type) tidal field £jk{t) inter- 
acts with the evolving quadrupole moment Ijk{t) of an isolated body, the tidal field does work on 
the body ("tidal work") — i.e., it transfers energy to the body — at a rate given by the same formula 
as in Newtonian theory: dW/dt = —^Ej^dlj^/dt. Thorne has posed the following question: In view 
of the fact that the gravitational interaction energy Ei^t between the tidal field and the body is am- 
biguous by an amount ~ Sjkljk, is the tidal work also ambiguous by this amount, and therefore is 
the formula dW/dt = —^£jkdljk/dt only valid unambiguously when integrated over timescales long 
compared to that for Ijk to change substantially? This paper completes a demonstration that the 
answer is no; dW/dt is not ambiguous in this way. More specifically, this paper shows that dW/dt 
is unambiguously given by ~^£jkdXjk/dt independently of one's choice of how to localize gravi- 
tational energy in general relativity. This is proved by explicitly computing dW/dt using various 
gravitational stress-energy pseudotensors (Einstein, Landau-Lifshitz, M0ller) as well as Bergmann's 
conserved quantities which generalize many of the pseudotensors to include an arbitrary function 
of position. A discussion is also given of the problem of formulating conservation laws in general 
relativity and the role played by the various pseudotensors. 
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I. INTRODUCTION AND SUMMARY 

For many gravitating two body systems in the uni- 
verse, the primary means of energy transfer from one 
body to the other is through tidal work. This work is 
accomphshed through the gravitational interaction be- 
tween the tidal field of one body and the mass multipole 
moments of the other body. A simple example of this 
is the work that the moon does on the earth as it raises 
the ocean's tides. Tidal work is also dramatically evi- 
dent in the moon lo, which gets heated as it travels in an 
elliptical orbit through Jupiter's tidal gravitational-field. 
This heating is the cause of lo's dramatic volcanism. In 
these cases it is clear that the tidal work is a physical 
observable and should in no way depend on one's means 
of calculating it. 

The term "tidal heating" is often used in place of 
"tidal work" , but is something of a misnomer. The net 
gravitational-energy that is transfered between two bod- 
ies interacting tidally does not necessarily go into heat. 
It may go into the energy needed to deform the body 
(i.e., raise a tide on it) or it may go into the internal vi- 
brational energy of the body. The net "tidal work" may 
also be negative, in which case the phrase "tidal cooling" 
might be more appropriate. Throughout this paper we 
will take the terms "tidal heating" and "tidal work" to be 
equivalent and to mean the net work done by an external 
tidal field on an isolated body. 

It seems evident that tidal work should be a "physical 
observable"; i.e., the net energy-transfer from one body 
to another should be a real, physical quantity and should 
not depend on the mathematics that one uses to calculate 
the work that is done. 

In calculating the tidal work for situations in general 
relativity, we consider the interaction between an isolated 
body and a complicated "external universe" (the precise 
definition of these terms will be discussed in Sec. III). In 
the most familiar cases, this external universe may simply 
refer to a gravitating body such as a star, a planet, or 
a black hole, that orbits around the isolated body. In 
such a system it has been shown by Thorne and Hartle 
that the total mass-energy M of the isolated body is 
ambiguous by an amount AM ~ XjkSjk, where X,fc is 
the mass quadrupole moment of the isolated body and 
£jk is the tidal field of the external universe. This mass 
ambiguity has the same physical origin as the ambiguity 
in the localization of energy in a gravitational wave: It 
arises from the fact that there is no preferred way to 
localize gravitational energy. This is true in Newtonian 
gravitational theory as well as in general relativity ||^ . 

This mass ambiguity shows up mathematically in the 
fact that the nonlinearity of Einstein's equations could be 
expected to produce a term ~ Ejkljk/i' in the time-time 
component goo of the spacetime metric outside the body 
(where r is a radial coordinate); and one is free, mathe- 
matically, to move this term or some arbitrary part of it 
into the mass M that appears in the standard equation 



500 = -l + 2Af/r + .... 

We can also understand this mass-energy ambiguity 
physically in terms of the standard experiment by which 
the total mass-energy M of a gravitating body is mea- 
sured: the application of the general relativistic ver- 
sion of Kepler's law to a test particle in orbit around 
the body. If the body is spherical and isolated and 
the orbit is circular, the body's mass-energy is related 
to the orbit's period T (as measured by distant clocks) 
and its radius r (defined to be its circumference divided 
by 2-k) by M = {r^ /G){2it /Tf . If the body is non- 
spherical, with various multipole moments including the 
quadrupole moment X^fe, then the moments perturb the 
orbit; but if one makes the orbit as circular as those per- 
turbations permit and measures the orbit's average ra- 
dius f, then the M that appears in the monopole 1/r 
part of the gravitational field is still given accurately, 
to first order in the moments, by the relativistic Kepler 
law M — (f^/G) (27r/r)^. The perturbations, being non- 
monopolar, disappear when averaged over the orbit. Sim- 
ilarly, if the body is precisely spherical but is perturbed 
by a weak, external tidal field £jk, then accurate to first 
order in those perturbations we can still compute M by 
this averaged-radius formula; again the perturbations av- 
erage to zero over the orbit (provided that the timescale 
on which £jk changes is long compared to the orbital pe- 
riod). However, if both a quadrupole moment Tjk and an 
external tidal field £jk are present simultaneously, then 
the product EjkXjk is monopolar in nature and has di- 
mensions of mass; and correspondingly Kepler's law with 
an average radius will give M + 0{£jk1jk)- Thus, one 
cannot measure M directly by Kepler's law. We shall 
discuss this mass ambiguity further in Sec. p^ . 

Zhang has used the Landau-Lifshitz pseudoten- 
sor (one of an infinite number of ways to localize 
gravitational-field energy) to derive the expression 
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for the rate at which a time-evolving tidal field £jk (t) does 
work on a body with time-evolving quadrupole moment 
Tjk(t). In view of the body's mass ambiguity AM ~ 
Ejkljk, Zhang (and also Thorne and Hartle ^) asserted 
that the work done should be ambiguous by an amount 
^ AM, and thus Eq. (|l|) should instead be written as 
dW/dt = {—^£jkdTjk/dt) and would be valid only when 
averaged over timescales long enough for W to build up 
by an amount large compared to AM = Ejkljk- This 
occurs, for example, in the long-term tidal heating of lo, 
during which £,jk and dXjk/dt oscillate partially in phase 
with each other, producing a cumulative work that goes 
into heat. 

More recently Thorne Q , while analyzing the effects of 
tidal forces on the stability of relativistic stars, claimed 
on physical grounds that Zhang |3| and Thorne and Har- 
tle [|| were wrong: The ambiguity AAf actually resides 
solely in the energy of gravitational interaction i?int be- 
tween the body and the external tidal field and not at 
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all in the body's self energy EscU (i-e. the total mass- 
energy contained within the volume of the body), and 
correspondingly not at all in the work done by the tidal 
field on the body, W = (change in -Egcif); so the rate of 
work done is unambiguously and instantaneously given 
by dW/dt = -\£jkdljkldt Q. 

An operational variant of Thome's argument is this: 
Consider a body on which tidal work is being done by the 
interaction between its time-changing quadrupole mo- 
ment 1jk{t) and some external time-changing tidal field 
£jk(t). One can imagine, at any moment of time, turn- 
ing off the tidal field while holding the body's size and 
shape unchanged (to first order in the tidal field). With 
the tidal field gone, we can imagine measuring the body's 
total mass-energy Mi, e.g., by the relativistic version of 
Kepler's laws. That measured mass-energy with tidal 
field (momentarily) gone can be regarded as the body's 
self energy itself (including rest mass). This measured self 
energy is unambiguous. Now, turn the tidal field back on 
and allow the system to evolve normally for some time 
At. Then, turn the tidal field off and make a second 
measurement of the body's total mass-energy Mf in the 
same manner as before. The difference Mf — Mi between 
these two measurements is the change in the body's self 
energy AEsdf- This is the work W done by the tidal 
field on the isolated body. We can then conclude that, 
since these measured changes in the body's self energy 
are unambiguous, dW/dt ~ dEscu/dt and W are also 
unambiguous. 

It is possible to test Thome's claim in a manner based 
on the following considerations: The self energy, defined 
in the above manner, will not change when the tidal field 
£jk{t) changes, but the shape and size of the body are 
held fixed and thence Ijk is held fixed. This is just 
a restatement of the fact that a force can do no work 
if there is no displacement. However, if Ijk changes, 
with £jk held fixed, then i?soif can change. This means 
that the unambiguous tidal work must be of the form 
dW/dt — dEsc\f/dt — (some constant) x £jkdljk/dt. 
The interaction energy, by contrast, should have the 
form of a product of the instantaneous tidal field and 
quadrupole moment, so its time rate of change should 
be a perfect time derivative dEint/dt — c?/dt[(some 
constant) x £jkljk]- The body's total mass M must be 
the sum of its self energy and that portion of the inter- 
action energy that resides inside and near the body (i.e., 
within the orbit of the test particle that one uses in ap- 
plying Kepler's third law to compute the mass), therefore 
dM/dt = dEscif/dt + dEi^t/dt ^ dW/dt + dEi^t/dt. If we 



express dM/dt in the form 

dKI . \ r> ctljk f ^ d\£']kX-']h 

—— = (constant) x £jk — ^ — h (constant) x — ^-r^^ 
dt dt dt 
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^Thorne Q needed this result as a key underpinning of his 
proof that tidal coupling stabilizes a star against gravitational 
collapse. 



then the first term must be dW/dt and the second 
dEint/dt. If Thorne is correct in his claim that dW/dt 
is unambiguous and that the total ambiguity of dM/dt 
resides in dEi^t/dt, then any computation of dM/dt us- 
ing any (general relativistically acceptable) localization 
of gravitational energy must give —1/2 unambiguously 
for the coefficient of the first term, while different local- 
izations should give different values for the coefficients of 
the second term. 

Purdue j|] has carried out detailed calculations of 
dM/dt in Newtonian theory using all possible localiza- 
tions of the gravitational energy and has found that, in- 
deed, the first term in Eq. (|^) always has the coefficient 
— 1/2 while the second depends on the localization. Pur- 
due has also verified that in general relativity, if one uses 
the energy localization embodied in the Landau-Lifshitz 
pseudotensor, but performs gauge transformations (in- 
finitesimal coordinate transformations) on the spacetime 
metric, the first coefficient (that associated with dW/dt) 
remains always —1/2, while the second (that associated 
with dEint/dt) changes with the changing gauge. 

In this paper we shall complete this test of Thome's 
claim: We shall verify that, when one changes the gen- 
eral relativistic energy localization by changing one's 
choice of pseudotensor, the first coefficient in Eq. (H) re- 
mains always —1/2 while the second changes and thus 
(partially) embodies the ambiguities present in localiz- 
ing gravitational-field energy. 

As a foundation for demonstrating this we first discuss 
m Sec. of this paper the problem of formulating co- 
variant conservation laws in general relativity. This is 
an underlying issue throughout this paper as the lack of 
an acceptable energy-momentum tensor for the gravita- 
tional field could possibly be a source of ambiguity in the 
calculation of the tidal work. We also discuss some of the 
various pseudotensors and conserved quantities that are 
used to describe gravitational-energy localization. 

In Sec. 01 we discuss the assumptions that go into 
our calculation of the tidal work. A key issue is that our 
calculations are performed in the local asymptotic rest 
frame, or LARF, of the body on which the work is being 
done. This means we are able to formulate only approx- 
imate conservation laws for our system. These laws are 
formulated in a buffer zone where the gravity of the iso- 
lated body is weak and the tidal field of the external 
universe is uniform. The spac etime metric of this buffer 
zone is described in Sec. Ill A. 



We then calculate t he tid al work using the Einstein 
pseudotensor ||^ (Sec. IIIB ) and review the calculation 
given by P urdu e H using t he La ndau-Lifshitz pseudoten- 
sor (Sec. [II C). In Sec. HID we perform the calcula- 
tion using the pseudotensor of M0ller which is sig- 
nificantly different from the two previously mentioned 
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pseudotensors, and in Sec. [HE we examine the calcula- 



tion using the conserved quantities found by Bergmann 
0. Bergmann's conserved quantities generalize many of 
the pseudotensors, including those of Landau and Lif- 
shitz, and Einstein. Each of these calculations gives the 
same, standard result dW/dt = —^£jkdTjk/dt for the 
tidal work, in agreement with Thome's assertion. 

Throughout this paper we adhere to the conventions 
of Misner, Thorne, and Wheeler |8j , hereafter referred to 
as MTW. Space-time indices are represented by Greek 
letters and spatial indices by Latin letters. We use units 
where G — c — 1. The constant on the right hand side 
of the Einstein field equations is -|-87r and the Minkowski 
flat-space metric rjap has signature (— , +, +, +). 



II. CONSERVATION LAWS AND 
PSEUDOTENSORS 

The formulation of covariant conservation laws has 
been a problematic issue since general relativity's for- 
mulation in 1916. The issue has been addressed by a 
large number of authors and some continue to work on 
this problem. 

If one considers a system without gravitational fields, 
as in special relativity, then the differential conservation 
laws for all matter and energy fields present are given by 
the familiar formula 



0, 



(3) 



where T^''' is the symmetric energy-momentum tensor of 
matter that appears as a source term on the right-hand- 
side of the Einstein field equations 



(4) 



By matter we mean all fields with the exception of the 
gravitational field. 

In general relativity Eq. (||) is not an acceptable con- 
servation law as it is not a tensor equation valid in all ref- 
erence frames. Instead we must use the covariant deriva- 
tive in place of the partial derivative and our equation 
becomes 



0, (5) 



where V^a-v are the connection coefhcients. From Eq. 

we can see that the mass-energy in matter fields is 
no longer conserved as energy can now be transfered be- 
tween the matter and the gravitational field. The quan- 
tity that is actually conserved in the sense of Eq. is 
some effective energy-momentum tensor T^^^ of matter 
plus gravitational fields which is given (in one variant) 
by Eq. (20.18) of MTW § as 



off 



(6) 



where t'^'^ is an energy-momentum pseudotensor for the 
gravitational field. In other variants, some of which are 



encountered below, Tl^^ = {-g)''/'^{Tt"' + t'"') where g = 
det\\gap\\ and n is a positive integer^. For each of these 
T^g , the equation T^^^-i, — can be rewritten as 



0, 



(7) 



and Tl^g can be written as the divergence of some "su- 
perpotential" H^^'^'^^ that is antisymmetric in its last two 
indices IsT: 



(8) 



Square brackets indicate antisymmetry of the tensor 
when the enclosed indices are swapped. Notice that Eq. 
(^) follows from Eq. (^) by differentiation and symmetry. 

As mentioned above, t^'' is not a true tensor, but rather 
is a pseudotensor that describes the localization of gravi- 
tational energy-momentum. That t^'' is not a tensor is a 
fact intimately linked with Einstein's Equivalence Prin- 
ciple. Since we are always free to choose our coordinates 
in spacetime to correspond to a freely falling frame where 
the acceleration vanishes at a point, we can equivalently 
choose a frame where the gravitational field vanishes at 
that point. In such a frame all the components of t'^'^ 
will likewise vanish at that point (provided one is using 
Minkowski coordinates). However, in any other reference 
frame, there is no reason why all the components of t^'' 
should vanish. Since any tensor that vanishes in one ref- 
erence frame must vanish in all reference frames, we can 
conclude that t^^" is not a tensor but a pseudotensor and 
quantities calculated from it will depend on the choice 
of one's coordinate system. To make matters worse, t^'^ 
is defined only up to a vanishing divergence, so there 
is an infinity of expressions for corresponding to an 
infinite number of ways in which one can localize the 
gravitational energy-momentum density. 

Despite their rather unpleasant nature in a theory so 
firmly rooted in the principle of general covariance, pseu- 
dotensors have proved to be rather valuable calculational 
tools, especially in gravitational-wave research (see for 
example [0). The reason is that, despite their nonco- 
variance, the T^^ can be used to compute covariant con- 
served quantities. For example, one can compute the to- 
tal 4-momentum of a system that resides alone in asymp- 
totically flat spacetime by the volume integral 



Ti:^d'x, 



(9) 



where d^x = dx^dx^dx^ is a 3- volume element of constant 
time. Even though the integrand depends highly on one's 
choice of coordinates, is a true vector that resides in 
the asymptotically flat region. 



■^We will sometimes refer to scalars, vectors and tensors with 
factors of {—g)"^"^ in front as scalar, vector and tensor densi- 
ties of weight n. 
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Using Gauss' law and the antisymmetry properties of 
the superpotential, it is also possible to express the 4- 
momentum as a surface integral: 



pp. 



rijd S, 



(10) 



where Uj is the unit normal vector to the surface S. It is 
important to note that these integrals must be evaluated 
using an asymptotically Lorentz coordinate system^ 



A. The Einstein Pseudotensor 

The first pseudotensor was formulated by Einstein in 
1916. The Einstein pseudotensor is often referred to as 
the "canonical" pseudotensor because it is derived using 
the general formula for the energy-momentum tensor of 
a classical field with Lagrangian density C and field vari- 
able r]A , which may be a tensor of any rank. In flat space- 
time this general formula is given by (see, e.g., Goldstein 
PI), 



dC 

driA,u 



C5, 



(11) 



and the Euler-Lagrange equations guarantee that 
T^^ = 0. In general relativity, the field variables are the 
components of the metric tensor y^j, , and the Lagrangian 
density is given by 



C 



Eq. (|llD then becomes 



f dC 



per p7 



0t) 



(12) 



(13) 



and again the Euler-Lagrange equations guarantee that 
{\/—gEtfi'^),i^ = in vacuum and e'J/j'^ 1^ ~ ^ where 



(14) 



when matter is present. The tensor density e2'^'' is often 
referred to as a "total energy-momentum complex" ; it is 
the Einstein variant of the T^^^ discussed above. 

From Eq. ( p^ we arrive at an explicit expression for 
the Einstein pseudotensor |13{: 



Note that raising or lowering an index for this pseudoten- 
sor does not produce a symmetric quantity, so we are 
unable to form a conserved angular-momentum complex 
from it. 

It was shown by von Freud @ that the Einstein com- 
plex can be written as the divergence of an antisymmetric 
"superpotential" pUL^'^^: 



(16) 



where 



{-.9(/V'-.9''V)},- (17) 



We can now form expressions for the covariant compo- 
nents of the 4-momentum of an isolated system by means 
of Eqs. (|) and (|l^: 



P„ 



(18) 



where the first integral is over the system's entire vol- 
ume, and the second is over a closed surface near spa- 
tial infinity (in the asymptotically fiat region of space- 
time). Because of the peculiarities of the pseudotensor, 
the above integral can only be interpreted as the covari- 
ant components of an energy-momentum 4-vector if one 
is using coordinates x'^ = {t,x,y, z) in which the metric 
5q/3 asymptotically approaches the Minkowski fiat metric 

To illustrate the coordinate dependent nature of the 
pseudotensors, we provide two well known examples 
(mentioned by M0ller and Anderson ) . If one were 
to calculate the integral / e^o*^ d'^x for the Lorentz metric 
5q/3 — Vaf3, its value would be zero, the expected energy 
for a region with no gravitational field. However, if we 
merely change to spherical coordinates, the value of this 
integral is infinite even though spacetime is fiat. Sim- 
ilarly, evaluation of the integral for the Schwarzschild 
metric only yields the mass M if one uses coordinates 
such that gap maps to rjap asymptotically as r ^ 00. 

Despite the restrictions on the use of this pseudotensor, 
it has still led to the reliable prediction by Einstein that 
gravitational waves exist and carry a definite energy. 



(15) 



^However, it should be noted that Nahmad-Achar and 
Schutz |ll| have devised a prescription for calculating 
pseudotensor-based conserved quantities for isolated systems 
in general relativity using coordinate systems with arbitrary 
asymptotic behavior. 



B. The Landau-Lifshitz Pseudotensor 

Landau and Lifshitz (LL) |^ were able to formulate a 
symmetric pseudotensor, thus allowing the construction 
of a conserved total angular-momentum complex. Their 
conserved total energy-momentum complex (their vari- 
ant of T^g ) is given by 



(- 



(19) 
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(where /I'^I'^'^l is defined below) and satisfies the usual 
property 



0. 



(20) 



The explicit form of t^L long and complicated. It is 
given by 

{~g)tZ^ = [Eq. (20.23) of MTW] (21) 

= [Eq. (96.9) of LL]. (22) 

The LL superpotential ft.^I'^'^l is related to that given 
by von Freud by : 



(23) 



The LL pseudotensor is related to the Einstein pseu- 
dotensor by the following formula ]l5[|: 



i-9Kl = {-9)9" V + iV^gn^^FU/^^l 



(24) 



spherical, but not to change the way one slices spacetime 
into space plus time. 

4. Under linear transformations, m'J^'' behaves like a 
mixed second-rank tensor. 

With the further restriction that it not contain higher 
than second order derivatives of the metric, M0ller explic- 
itly exhibits a unique energy-momentum complex with 
these properties, and he shows that it can be written as 
the divergence of the following antisymmetric superpo- 



tential 



AM 



dx 



1 



d 



8tt dx" 



^ig^c.,0-9^.0,o.)g''^9^"]■ (27) 



The M0ller superpotential 
Freud superpotential by |15|, 



l/cr] 



is related to the von 



X, 



^FU- ' (28) 



As in the case with the Einstein pseudotensor, integrals 
of the Landau-Lifshitz pseudotensor also produce strange 
results in curvilinear coordinate systems. Asymptotically 
Lorentz coordinates must again be used if one wants sen- 
sible results. 



C. The ]V[0ller Pseudotensor 

The M0ller pseudotensor is significantly different 
from the two complexes mentioned above. In deriving 
his pseudotensor M0ller sought to eliminate the problem 
that the integral given in Eq. (|lj) yields strange results 
if one converts to curvilinear coordinate systems. 

To define his pseudotensor, M0ller makes use of the 
fact that one can always add a quantity S'p'^ to the Ein- 
stein complex and still retain energy-momentum conser- 
vation, provided that 5^"^ ^ = 0. The new total pseu- 
dotensor complex (matter plus gravitational fields) will 
thus have the form 



M- 



(25) 



M0ller additionally restricts the form of m^^^"^ by requir- 
ing the following conditions (see also Komar p6[): 



1. It must be identically conserved: yi'^fj 



= 0. 



2. The integral over some 3- volume of constant time 
of u'^ti" must produce the same results (in an asymp- 
totically Lorentz coordinate system) as Eq. (Eq): 



M 



(26) 



Like the Einstein pseudotensor, the M0ller complex is 
not symmetric and thus cannot be used to form conserva- 
tion laws for angular momentum. Moreover, unlike the 
complexes of Einstein and Landau-Lifshitz, the M0ller 
complex is not entirely quadratic in the first derivatives 
of the metric but has a term that is linear in the second 
derivatives of the metric. As pointed out by M0ller him- 
self Q, this means that M^f/ will generally not vanish 



in a local Lorentz frame with no matter present 
will become an issue in Sec 



[1] 



This 

when we use the M0ller 



complex to calculate the tidal work. 



D. The Bergmann Conserved Quantities 

Recognizing that conservation laws are related to the 
invariance properties of physical laws, and combining 
this with the fact that the equations of general relativity 
are invariant under arbitrary coordinate transformations, 
Bergmann proposed that to each infinitesimal coordi- 
nate transformation there would correspond a conserved 
quantity. Making use of various identities, Bergmann 
[|7| constructs a relationship between an arbitrary vector 
field (which may be thought of as inducing infinites- 
imal coordinate transformations on the metric) and the 
generators C of these transformations: 



^gG^^Sg, 



where 



(29) 



(30) 



3. m'^o^ and m'^o'^ behave like scalar and vector den- 
sities under arbitrary changes of the spatial coordinates. 



^new 



^ \-^old ' ■''old ' -^old A -^ne 



■-old- 



This allows one 



to change the coordinate system from say, Minkowski to 



Bergmann chooses C = 2-^— gC'^^o- as a solution to 
Eq. (E9). However, one may always add an arbitrary curl 
field V^P''\„ to CP and still satisfy Eq. dggj). Bergmann 
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chooses this curl such that the resulting expression con- 
tains no higher than first derivatives of the metric. His 
final expression satisfying Eq. (|9|) is 

c'' = 2r\/^G.^ + (rFC/[^"i) . (31) 

Eq. ( |3l| ) represents a weakly conserved quantity, mean- 
ing that it satisfies C'^ = whenever the vacuum field 
equations are satisfied (C^"^ = 0). The corresponding 
strong conservation law is D^^ — 0, where 

^ ^'^ (32) 

as can be easily shown if we make use of the antisym- 
metry of T^U^^^ and the commutativity of partial deriva- 
tives. 

From this strongly conserved quantity Bergmann con- 
structs the 4-momentum in the same manner as we did 
fromEq. (|): 

= J D^'d^x. (33) 

Using Gauss' law and the antisymmetry of pJ/i^''', we 
can write the fi — component as a surface integral. 




According to Bergmann, from the weakly conserved 
quantities C^, expressions equivalent to several of the 
pseudotensors can be derived by making specific choices 
for For example choosing = k"' yields the 

canonical-Einstein expression contracted with k"^ , while 
setting = \/~gg'^°'ka yields the Landau-Lifshitz ex- 
pression contracted with k'^ . We have been unable to 
find a similar choice that yields the M0ller pseudotensor. 
We believe that this is due to Bergmann's choice of the 
curl field V^P'^\a- containing no second derivatives of the 
metric. Recall that unlike most other pseudotensors, the 
M0ller complex contains second derivatives of the metric. 

E. Other Formulations of the Conservation Laws 

Aside from the method of using pseudotensors to for- 
mulate conservation laws in general relativity, there exist 
several other approaches as well. One of these is the 
method of quasilocal energy, a covariant definition of en- 
ergy that arises from a Hamiltonian formulation of gen- 
eral relativity. For a discussion of the equivalence of the 
quasilocal and pseudotensor approaches to gravitational 
energy-momentum, see Chang, Nester, and Chen | [l7| . In 
parallel with this present research. Booth and Creighton 
have calculated the tidal work using the Brown and 
York |l9|| quasilocal energy approach, and have arrived at 



the same result as we deduce below using various pseu- 
dotensors: dW/dt = —^£jkdJjk/dt. 

There have also been efforts to find an energy- 
momentum tensor for the gravitational field. Shortly 
after relativity was formulated, Lorentz and Levi-Civita 
proposed that the Einstein tensor G^" be used as a gravi- 
tational energy-momentum tensor. This however, did not 
prove fruitfulQ Most recently, Babak and Grishchuk ]2l| ] 
have shown that when one formulates general relativity 
as a non-linear field theory in fiat spacetime, then there 
exists an energy-momentum tensor for the gravitational 
field that has all the nice properties one might wish. This 
energy-momentum tensor is the Landau-Lifshitz pseu- 
dotensor with partial derivatives replaced by covariant 
derivatives with respect to the flat-background metric. 

III. CALCULATION OF THE TIDAL WORK 

In calculating the tidal work, we consider a system 
consisting of an isolated body that interacts with a com- 
plicated external universe in the slow-motion approxima- 
tion. The body is isolated in the sense that the radius of 
curvature TZ of the external universe and the lengthscale 
C on which this curvature changes must both be large 
when compared with the size R of the isolated body: 
R/TZ ^ 1 and R/C <C 1. This means that the exter- 
nal universe is not subjecting the isolated body to very 
strong gravitational fields (as would happen e.g., in a 
neutron star and black hole close to merger) and that 
the tidal field of the external universe is nearly uniform 
in the region near the isolated body. By slow-motion, we 
mean that the timescale r for changes in the mass and 
current moments of the body and the tidal field of the 
external universe are small compared to the size of the 
body: R/t ^ 1. If this were not the case, we would 
have to worry about changes in the mass-energy M due 
to gravitational radiation and other higher-order effects. 
For detailed discussions of the constraints R/TZ ^ 1, 
R/C <C 1, and R/t <C 1, and of various approximations 
based on them (which we shall use below), see Thorne 
and Hartle and the recent paper by Purdue , whose 
analysis we are continuing. 

Some examples of isolated, slow-motion bodies dis- 
cussed by Purdue ||^ include: (i) a compact object such 
as a neutron star or black hole in a binary inspiraling 
system that is not too close to merger; and (ii) Jupiter's 
moon lo, which gets tidally heated as it travels through 
Jupiter's tidal field in an elliptical orbit. 



*For an excellent discussion of the exchange that occurred 
between Einstein, Levi-Civita, Lorentz, and others concerning 
conservation laws and the prediction of gravitational waves, 
see the article by Cattani and De Maria pOjl 
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Our calculation of the tidal work involves computing 
dM/dt, the rate of change of the mass of an isolated body, 
and then expressing dM/dt in the form of Eq. ^ and 
reading off the two coefficients. We use the multipole mo- 
ment formalism discussed in Thorne [ p2[ and Thornc and 
Hartle and treat gravity as a non-linear field theory in 
flat spacetime. The computation of dM/dt is carried out 
as a 2-dimensional surface integral of a pseudotensor in 
the "buffer zone" or local asymptotic rest frame (LARF) 
of the isolated body (see Eq. (2.3a) of Thorne and Hartle 



dM 
~dt 



^^^it^'d^S 



■'0 1 



(35) 



where d^Sj = rijr^dil is the surface element of a 2-sphere 
dV in the buffer zone that encloses a volume V and has 
unit normal rij and solid angle dfl. This buffer zone is a 
region that surrounds the isolated body but is far enough 
away that gravity in it can be considered weak. At the 
same time, it is close enough to the body that the tidal 
field of the external universe appears homogeneous. The 
buffer zone can be described more precisely 1^ as the 
region where r/C <C 1, r/TZ <^ 1, and M/r <C 1, r being 
the radial distance from the isolated body. The rate of 
change of mass-energy through the surface dV is dM/dt, 
and M is the total mass-energy inside V. Note that our 
analysis is thus valid even for a strongly gravitating body 
such as a black hole, provided there exists a buffer region 
around it where gravity is weak, the external curvature 
is nearly uniform, and the spacetime curvature is not 
changing too rapidly. 

For the purposes of our discussion and to the order that 
our calculations are valid, there are only three relevant 
parameters that characterize the spacetime: 

1 . The total mass-energy M of the isolated body. 

2. The quadrupole moment 2jk of the isolated body, 
which, in the limit of weak gravity, is given by 



Ljk 



1 



pXjXk - ^SjkT^ d^x 



(36) 



3. The tidal field Sjk = Rjoko of the external universe, 
where Rap-yS is the Riemann tensor of the external uni- 
verse. 

Note that both Xj^. and £jk are symmetric and trace 
free tensors that reside in the buffer zone, and that we are 
using coordinates that are as Lorentz as possible (with re- 
spect to the physical metric) throughout the buffer zone; 
i.e., 5a/3 = Vc.f3 + 0iM/r) + 0{I/r^) + 0i£r^). These I^k 
and £jk are spatially constant in the buffer region but 
they may depend on time. 

The body has additional multipole moments: the cur- 
rent quadrupole moment Sjk, the mass octupole moment 
Tjki, etc; and the external universe has additional tidal 
fields of "magnetic- type" {Bjk, ■ ■ ■ ) and "electric-type" 
{£jkh ■ ■ ■ )j see e.g., Thorne and Hartle [Q. These mo- 
ments and tidal fields can couple to each other to produce 
tidal work: dW/dt ~ BjkdSjk/dt Sz Sjkidlju / dt & • • • . 



In some situations these contributions to dW/dt might 
be larger than the one, dW/dt £jkdXjk/dt, that we are 
studying, but typically the mass quadrupole will dom- 
inate. In this paper we restrict ourselves to the mass 
quadrupole term, whether or not it dominates the tidal 
work, because we are seeking to discuss an issue of prin- 
ciple first raised by Thorne and Hartle ||l|: the non- 
ambiguity of the £jkdXjk/dt tidal work. Presumably our 
results can be generalized to higher order moments, but 
we shall not attempt to do so. Correspondingly, in our 
analysis we shall consider only M, Tji^, and £jk- 

Keep in mind that when we identify a mass M as in 
Eq. (35), we are only doing so in an approximate sense. 



This is because the mass, momentum, and angular mo- 
mentum only have precisely defined values in an asymp- 
totically flat spacetime. As our spacetime is only locally 
asymptotically flat, the conservation laws only give ap- 
proximate values of the mass, momentum, and angular 
momentum in the buffer zone where spacetime is approx- 
imately asymptotically flat. 

In particular, as we discussed in Sec. ||, there is an am- 
biguity in the mass M of the isolated body AM ^ Tjk£jk- 
If Xjk is oscillating in time, then this ambiguity is of the 
same order as the amount of energy that is transfered 
between the isolated body and the external universe by 
tidal work during one period of oscillation. We can un- 
derstand this ambiguity more clearly by examining the 
time-time part of the metric in the buffer zone IS]: 



2M 

goo = - 1 H 1- Slab 

r r 



x'^x'' 



(37) 



where r is the distance from the center of the isolated 
body as measured in its local asymptotic rest frame. We 
have omitted terms involving higher order mass and tidal 
multipole moments (e.g. octupole moments Xj^i and 
£jki) and also terms that are products of M, Xjk, and 
£jk which result from the nonlinearities of the Einstein 
fleld equations. One of these nonlinear terms has the 
form 



(38) 



which has the same form as 2M/r, i.e. (monopole)/r, 
and which has a coefficient that is gauge dependent. The 
similarity in structure between M/r and (Xjk£jk)/r im- 
plies that it is possible to move portions of the gauge- 
dependent term given in Eq. ( p^ ) into or out of the 
2M/r term. One can interpret this as meaning that the 
mass AI that one reads off the metric is ambiguous by an 
amount on the order of 



AM - X,k£ 



(39) 



Purdue Ql shows that this ambiguity is also present in 
Newtonian theory in the form of an ambiguous gravita- 
tional interaction-energy inside and near the body. More 
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specifically: The total mass-energy Ey enclosed in the 
volume V can be expressed in Newtonian theory as 



Ev — E, 



self - 



E, + Ei 



int 5 



(40) 



where Eseii is the isolated body's self-energy (which de- 
pends on the body's rest-mass and internal energy den- 
sity distributions), E^. is the external field energy inside 
the volume V (which depends only on the external tidal 
field £jk ) , and Emt is the interaction energy inside V and 
is given by 



Eint — 
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(41) 



Here the coefficient a depends on one's choice of Newto- 
nian energy- localization. For example (see Purdue [^), 
the choice a = localizes all of the gravitational en- 
ergy in the field, so that the total gravitational energy is 
given by a volume integral of (V$)^/(87r) (which should 
be familiar from electrostatics) . Alternatively, the choice 
a — 1/2 localizes the gravitational energy entirely in the 
matter, so the total-gravitational energy is given by a 
volume integral of Here p is the mass density, so 

clearly it vanishes outside the material of a gravitating 
body. 

From Eq. (^) we can see that a Newtonian calculation 
of dEv/dt (which is analogous to our general-relativistic 
calculation of dM/dt) will include a term that is the 
total time derivative of the interaction term Ei^t- Our 
general-relativistic calculation will also have an interac- 
tion term of this same form, where the coefficient in front 
depends on the energy-localization scheme. In the New- 
tonian case, Purdue shows that, despite the ambiguity 
of dEint/dt, the rate of change of the body's self energy 
dEscif/dt is given unambiguously by the tidal work for- 
mula dEsou/dt = —-^SjkdXjk/dt. Our general relativistic 
analysis will produce this same conclusion. 



A. Metric in the Buffer Zone 

In our calculation we will consider general relativis- 
tic gravity not as a geometric phenomenon involving the 
curvature of spacetime, but rather as a non-linear field 
theory in flat spacetime. We treat the field variables 
(the metric components gap) as a perturbative expan- 
sion in some dimensionless parameter e which is actually 
the gravitation constant, G = 1 in our system of units. 
Thus, terms of 0{e) are linear perturbations around flat 
spacetime; terms of 0{e^) are quadratic, etc. All raising 
and lowering of indices is done with the flat Minkowski 
metric r]af3. 

Our three parameters M, Tjk, and £jk can all be con- 
sidered linear in e. In our calculation of the tidal work, 
it is clear from the form of Eq. (^ that we may also 
need to consider terms in the metric that are quadratic 



in e, as such terms may go as TjkSjk- We thus expand 
the metric up to quadratic order: 



gap = Vap + Shap + E^kap, 



(42) 



where hap contains terms that are linear in M, Ijk, and 
Ejk and their time derivatives, while kap contains terms 
that are products of any two of those three quantities 
and their first time derivatives (for example MX, XX, 
Xf , ££, etc). Any terms that are cubic or higher in the 
perturbation expansion e cannot contribute to the tidal 
work (which is itself of order e^) and can be discarded at 
any point in the calculation. 

Since the mass quadrupole moment Ijk and the tidal 
field £jk are spatially uniform in the buffer zone, the spa- 
tial gradients of these functions vanish. Furthermore, 
since dt ~ 1/t, the slow-motion approximation allows us 
to ignore all second and higher order time derivatives of 
Xjfc and Ejk- 

We impose separate gauge conditions on the linear hap 
and quadratic kap parts of the metric. The linear part 
of the metric has been calculated by Zhang a in the de 



Bonder gauge and is used by Purdue | 
The de Bonder gauge in linear order 



I in her analysis, 
is given by the 



condition that h^^ ^ - 
metric perturbation: 



0, where h^^ is the trace reversed 



2 ' 



(43) 



In terms of the metric perturbation hap we can write this 
as the condition: 



where h is the trace of hap'- 

h = ha''= r^hap. 



(44) 



(45) 



To the order we should need in the slow-motion ap- 
proximation (i.e. neglecting second and higher order time 
derivatives) and ignoring higher order multipolc moments 
(i.e. Tjki, Sjk, Sjkii Bjk, etc.) the linear part of the 
metric is given in Cartesian coordinates and de Bonder 
gauge, accurate to 0{e) in the buffer zone, by 0: 



hoo = -2$ - 2iji + 3^:^^ - 

ijaX^ 10 • a b ' 4 ■ a 



r3 21 



21 



h,, = -2$5,, = 6,,i2— + S^jx'x^), 



(46) 

2 

(47) 
(48) 



The general de Bonder condition is given by g*^ 



where g*^ 



-99'^ 
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where $ is a scalar potential analogous to that of New- 
tonian gravity and Aj is a vector potential that has no 
Newtonian analog. 

This metric, by virtue of the dc Bonder gauge condi- 
tions and the approximations mentioned above, satisfies 
the following relations, which wc shall use in our calcu- 
lations below: 



^i.o = 0, 



= -4$,o, 
0. 



ha"" 



(49) 
(50) 
(51) 
(52) 



To compute the quadratic part of the metric kap one 
would have to solve the second order Einstein field equa- 
tions. It turns out that for the pseudotensors we are 
considering, a suitable choice of gauge will make the di- 
rect calculation of kap unnecessary. This specific gauge 
will be discussed later. 

Other formulas that will be useful in the calculations 
that follow are given below, accurate up to order e^: 



(53) 



(54) 



- 5 = det(5a/3) ^l + eh+ -e^^h^ - h'^'^Kp) + e^k. 

2 (55) 

Taylor expanding Eq. ( p5[ ) about e, we also have, accu- 
rate to order e^, 



(56) 



(57) 



We are now ready to compute the tidal work using vari- 
ous pseudotensors. 



B. Calculation of Tidal Work using the Einstein 
Pseudotensor 

We wish to calculate the integral given by Eq. (|3^) for 
the rate of change of mass-energy of the isolated body 
using the Einstein pseudotensor^— g^t^^ : 



dM 
~dt 



-gEto^rijr'^dfl, 



(58) 



where Uj = Xj/r is the unit normal to a surface lying 
in the buffer zone at some radius r, and dO. is the solid 



angle on that surface. If we look at the form of \/—gEti_i'^ 
given by Eq. (15), we can see that in order to obtain an 
expression that is accurate up to order we only need 
to expand the metric to order e. This means that we 
can ignore all terms appearing in Eqs. (^), ( [s^ ) - ( [S^ ) 
that go as e^. Eq. (15) thus takes the form (in a general 
gauge): 



(59) 



1 



'-h°"'''hc 



1 



Applying the de Donder gauge constraints (^), we see 
that the first and second and the seventh and eighth 
terms cancel each other. 

From Eq. (|58| ) we see that we only need to evaluate 
the E^is' terms. These terms evaluate]^ to 



-9eW 



1 

47r 



(60) 



Using Eq. (46) we see that 



15 • X qXi^x qX dX j 

~7~-^ab^cd 



r' 

Xq^XJjXq 



3^ 
2- 



X(iX}jX(2 



(61) 



Note that we have ignored terms that go like MM, MX, 
M£, TM, XT, £M, and ££ as they do not contribute to 
the tidal work. That these terms do not contribute is ap- 
parent if one considers that the tidal work must arise due 
to a coupling between the mass multipole moments of the 
isolated body and the tidal field of the external universe. 
This, combined with dimensional considerations, implies 
that only terms of the form 2£ and X£ can contribute to 
dW/dt. Also we do not worry about spatial indices being 
up or down since we are using Cartesian coordinates. 

We must now evaluate the surface integral given in 
Eq. (M). Note that keeping our calculations accurate to 
order justifies our setting the factor of on the left- 
hand-side of Eq. (^) to unity. To perform the surface 
integrals over the terms in Eq. (61), we first note that 
since the multipole moments do not vary spatially in the 
buffer region, they can be pulled out of the integrals. The 
surface integrals that remain are all of the form 



nanbUc ■ ■ - Up dfl , 



(62) 



''Part of this calculation was performed using the tensor al- 
gebra package GRTensorll ]23|. 
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where = x^/r is a component of a unit radial vector. 
Evaluating such integrals (see Sec. IIB of Thorne 22|) 
we finally arrive at 



dt [ lO^^'^^^" 



(63) 



We can identify the second term in the equation above 
as the analog of the Newtonian interaction-energy given 
by Eq. (^ ) where the Einstein pscudotensor localization 
corresponds to the choice a — 1. That this term is in fact 
the derivative of an interaction-energy term is apparent 
if once considers that the interaction energy Ei^t between 
an isolated body and the external universe must depend 
only on the instantaneous fields and can only be given by 
a product ~ Tjk£jk ■ The rate of change of this interaction 
energy must then be a perfect differential. Also, a term 
that goes like ~ JjkSjk could not contribute to the tidal 
heating since no work is done if the isolated body does not 
change (just as no work is done if a force is exerted but no 
displacement results); see Sec. I for further discussion. 
These facts indicate that the first term in Eq. ( |63| ) is the 
tidal work while the second term is the rate of change of 
the interaction energy between the external universe and 
the body; cf. Eq. (||). 



C. Calculation of Tidal Work using the 
Landau-Lifshitz Pseudotensor 

The calculation of the tidal work using the Landau- 
Lifshitz pseudotensor is very similar to that for the Ein- 
stein pseudotensor shown above. This calculation was 
performed by Purdue ||^ and we will only summarize her 
results here. 

Since the Landau-Lifshitz pseudotensor, like the Ein- 
stein pseudotensor, is quadratic in the first derivatives of 
the metric, one only needs to expand gafj to first order 
in e. We thus only need to consider the linear part of 
the metric ha/3 when evaluating the integral in Eq. ( |5^ ) 
(where ^/^stfj.'^ is replaced by {—g)t'^'i)- 

Evaluating this integral and keeping only terms that 
contribute to the tidal work (in the same manner as the 
previous section), Purdue arrives at 



dt ^ 2 ^'^ dt 



d 
dt 



-It f 



(64) 



We notice that changing the energy localization scheme 
from that of Einstein to that of Landau-Lifshitz has sim- 
ply changed the coefficient of the second term, which 
we again identify as the derivative of the interaction en- 
ergy. Note that the Landau-Lifshitz localization scheme 
is analogous to a choice of a = —3 in the Newtonian in- 
teraction term. We also see that the tidal work term (the 
first term in Eq. (|6^)) has remained unaffected. 



D. Calculation of Tidal Work using the M0ller 
Pseudotensor 

We shall now perform the calculation of dM/dt once 
again, this time making use of the M0ller pseudoten- 
sor. Since we are working in the vacuum buffer zone 
where T'^'^ = 0, we can use Eq. (27) as the expression 
for the energy-momentum pseudotensor of the gravita- 
tional field. Note that Eq. (27) is actually the total con- 
served complex that we would use if non-gravitational 
fields were also present. 

If we examine closely the form of Eq. (27) we will find 
that unlike the two previous pseudotensors discussed, the 
M0ller complex has a term that is linear in the second 
derivatives of the metric perturbations. This means that 
we will not only have terms like h"^ ^f^hap''^ but will also 



have terms like hn 



and fc^'^ 



This means that it is 



important that we expand the metric up to quadratic 
order in e as these terms that are linear in the sec- 
ond derivatives of the metric perturbation are actually 
quadratic in s and will thus contribute to the calculation 
of the tidal work. 

Using the metric given by Eq. (^) in Eq. (27), we 
arrive at the form of the M0ller total energy-momentum 
pseudotensor, correct up to order e^, in a general gauge: 

+ \e\Kh^-'- + hh^'^/ - h^h/-'^ - hh^\f) 

(65) 

If we work in de Donder gauge to linear order we can use 
Eqs. (46) - (48) for h^p ■ However, we still do not know 
the form of kap. 

Fortunately, we can make use of the vacuum Einstein 
field equations 

Rfzi/ — r fzi/^a r ^a^y ~\~ L (Pa^^ E ^cl — 

(66) 

to solve for the derivatives of kap as they appear in Eq. 
(65). Specifically: If we substitute the metric ( ^ ) in Eq. 
(|66|) we can expand the Ricci tensor in powers of e: 



[h] + £2i?(2j [hh] + e2i?(2j [fc] ^ 0, 



(67) 



where the superscript on the W^J means that the in- 
dicated piece of the Ricci tensor contains only terms of 
order e". The terms in the brackets indicate that the 
part of the Ricci tensor in question contains terms that 
go like the indicated multiple of the metric piece {h^p or 
kap) and its derivatives. Eqs. ( [70| ) and (71) make this 
clear. We now require that the vacuum field equations 
vanish in each order of e: 



i?«[/i]=0. 



(68) 
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(69) 



The first of these equations yields the linearized vacuum 
field equations: 



h,fj.i, = 0, 



(70) 



(2,) 

while in the second (|69|) one can solve for RliJlkap] to 
give: 



ua 4_ J." _ J, 



(71) 

We can now use Eq. (71) to substitute for the kafs terms 
that appear in Eq. (65), if we pick a gauge for the kajs 
such that 



1^ fi,ua ^ i^,/^a 



This is easily done if we choose the gauge 



k " — h 



(72) 



(73) 



We can now use Eqs. (71) and ( Jt^ ) along with the 
de Bonder gauge conditions and Eqs. (51) and (52) to 
simplify the M0ller pseudotensor (65) into the form: 



1 



uaa,vi _ uaa, i- 



Inserting Eqs. (46) - (48) for hap, we obtain!: 



(74) 



(75) 



Expanding this expression and ignoring terms in the 
same manner as Eq. (61), we plug into the integral 



dM 
~dt 



M 



and arrive at 



dM . . 



(76) 



(77) 



where the M term comes from the — 2$jo term in (75), 
which in turn arises from the fact that the M0ller pseu- 
dotensor contains a piece linear in the second derivatives 
of the metric. Eq. ( [77| ) can then be written as 



dt ~ 2 dt ' 



(78) 



Note that in the case of the M0ller pseudotensor, the 
perfect differential term that represents the interaction 
energy vanishes. This serves to support our intuition that 
the interaction term is simply a mathematical artifact of 
our choice of energy localization and that the tidal work 
is in fact uniquely given by Eq. (|^). Comparison with Eq. 
( ^ ) implies that the M0ller pseudotensor corresponds to 
a Newtonian energy localization given by a = — 2. 



E. Calculation of Tidal Work using the Bergmann 
Conserved Quantities 

To further confirm our intuition about the tidal work, 
we can calculate dM/dt using the Bergmann conserved 
quantities given by Eq. (^2[). However, before we do this, 
we must determine the form of the arbitrary vector field 
on which the momentum density D^^ depends. This 
form must be such that the integral (|3^ ) gives 



eM + 0{s^), 



(79) 



the mass-energy M plus terms of order which arise 
from the fact that the surface of integration dV lies in 
the buffer zone, where spacetime is only locally asymp- 
totically flat. 

Computing via Eq. (^J) first involves calculating 
the von Freud superpotential (17) in terms of the metric 
expansion given by Eq. (^2|). After a rather laborious 
computation, we find the von Freud superpotential in 
general gauge and to quadratic order in the metric per- 
turbation: 



-h h. 



+ {5cPkP\x - 5jk'^\x) + {kc?'^ - fc/'^) |. (80) 

As in the case of the M0ller pseudotensor, we have 
terms that are linear in the derivatives of the kafi- For 
our calculation of dM/dt below, we will find that these 
terms may again be expressed in terms of the hap by the 
same choice of gauge as we used in the previous section, 
Eq. (|73|). However, we cannot do this for our present 
calculation of P'^ . Fortunately, this will not be a problem 
since we only wish to show that reduces to the mass- 
energy M plus terms of order Since kap is of order 
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it is not necessary to include the terms containing kap in 
our calculation of the surface integral in Eq. . This is 
also true of the terms in (80) that are products of ha/3 and 
its derivatives. We therefore only need concern ourselves 
with the first six terms in Eq. (80) (those linear in haf3 
and its derivatives) as only these terms could conceivably 
affect P° at order e. 

Now, let us make a guess as to the form of the vec- 
tor field on which depends. A volume integral 
of Bergmann's (or the equivalent surface integral of 
(^)) will reduce to Eq. (|7^) only if the form of ^'^ is prop- 
erly constrained. If our spacetime were precisely asymp- 
totically flat, we would expect to be asymptotically 
the timelike Killing vector field, d/dt + 0{l/r). 

Since our spacetime is not asymptotically flat but only 
locally asymptotically flat, it is reasonable to expect that 
is locally asymptotically Killing, by which we mean 
that it can be written as the sum of a timelike vector, 
6q , plus deviations C"^, from this timelike vector that are 
due to the fact that the spacetime is not flat: 



(81) 



As for the form of C^, we can construct a quantity 
which is the most general vector field that is a) dimen- 
sionally correct and b) constructed only out of the param- 
eters that characterize our spacetime in the slow-motion 
approximation: M, X,fe, £jk and their first time deriva- 
tives. Such a vector field has the following form (to order 
e): 



( = aM + b h C 



r 

j-^ab-^a-^b 



^^ab-^a-^b ~t~ f^ab'^a-^b'^^ 

(82) 



Mxk , ,Mxk , labXaXbXk , .labXaXbXk 

^ = 9—;^ + + 1 +J 

, 7 -^akXa , j-^akXa , ^abXaXbXk 

+ k J 1- i g h m 

-I- nSabXaXbXk + oSakXar + pSakXaT^ , (83) 

where the coefficients a through p are real-valued con- 
stants. 

Using this prescription, we are now able to evaluate P^. 
Our result reduces to Eq. (|79|) as required. We thus see 
that the surface integral (piFdoes indeed give the mass 
as the only 0(e) term, plus terms of O(e^) that arise from 
the facts that our spacetime is locally asymptotically flat 
and that we evaluate P° on the 2-surface dV that lies at 
some finite r in the buffer zone. 

Now we wish to compute the tidal work, dP'^/dt or 
dM/dt. Taking a time derivative of Eq. (|33| ) and apply- 
ing Gauss' law, we easily arrive at the expression 



(84) 



Plugging Eq. (^l|) into Eq. (^) allows us to write 



dM 

IT 



FUQ''\^d'^Sj - 



(85) 



Let us now think carefully about what terms we actu- 
ally need to calculate in these integrals. If we examine 
the first integral in Eq. ( p5[ ) we realize (comparing with 
Eq. (80)) that all of the terms in the integrand are of 
order and can therefore contribute to the tidal work. 
However, we again have the problem that we do not know 
the explicit form of kaf}. Fortunately, if we expand the 

terms in (fJ/ct"'') ^ that depend on kajs and apply the 
gauge condition of Eq. (fz^), we find that these terms 
reduce to 



a .y (7,1/ } ' 



(86) 



Comparison with Eq. ( fz^ ) shows that this is the same 
situation encountered with the M0ller pscudotensor. We 
can thus use Eq. (71) to express the first integrand in 
Eq. (^) entirely in terms of hafi and its derivatives, just 
as we did in the previous section. Applying de Bonder 
gauge and a fair amount of algebraic manipulation we 
arrive at 

- hxp.^hf^P^^ + h^P^f^hxp.,, - S^'^h^'^'-'h^p... (87) 

Evaluating the comnonents we need for the first integral 
in Eq. (p5|), we findEl 



If we expand this expression as we did in Eq. (61), keep- 
ing only terms that can possibly contribute to the tidal 
work (the products of 2£ and If), and evaluate the first 
integral in Eq. (^), we have 



dljk 
~dt 



--{I,k£,k) 
10 dt ^ 



(89) 



One should note that this is precisely what we obtained 
using the Einstein pscudotensor (Eq. (63)), and indeed, 
the left-hand side of (|9|) is the surface integral that one 
would obtain by substituting Eq. (16) into Eq. (^ 

We must now evaluate the second integral in Eq. (|8£ 
to see if it will contribute to the tidal work. The compu- 
tation of this integral is made much simpler if we realize 
that (^^ is linear in e. Since the tidal work is quadratic 
in e we only have to concern ourselves with the piece of 
pUa'^^ that is linear in e. Applying de Bonder gauge 
to the linear piece of Eq. (80) we can write the second 
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integral in Eq. as 



1 

1^ 



+ 0{e'), (90) 



where we have apphed Eqs. (51) and (52). A few pages of 
algebra, and again ignoring terms that cannot contribute 
to the tidal work, reduces this to 



/(CpC/ri)^.'5,.(Ao+2.)|,X,..5,.) 



(91) 



Combining Eqs. ( |89| ) and (^) we finally arrive at an 
expression for dM/dt in terms of the products that can 
contribute to the tidal work (I£ and T£): 



dM 
IT 



1 dijk 
'2^'" dt 



3 2 3 \ d , ^ , 

10 + -5'+ rjdt'^^^'^^'l^^^ 



We can now see conclusively that changing the energy 
localization has no effect on the tidal work term and 
merely changes the coefficient in front of the arbitrary 
interaction-energy term. 



IV. CONCLUSIONS 

This paper completes a demonstration that the tidal 
work caused by the interaction of an isolated body's 
quadrupole moment Ijk with the electric-type tidal field 
£jk of an external universe is unambiguous, despite the 
ambiguity in the definition of the mass M of such a sys- 
tem. Purdue demonstrated that a gauge change does 
not lead to any ambiguity in the tidal work in general 
relativity, nor does one's choice of energy localization in 
Newtonian theory. We have shown that this energy local- 
ization invariance carries over into the general relativistic 
description. In addition, the work of Booth & Creighton 
Ip^ , carried out simultaneous with our own, supports the 
conclusions of Purdue and ourselves through an indepen- 
dent, though equivalent, approach using the techniques 
of quasilocal energy. 

One of the main motivations for this paper was to 
strengthen the arguments used by Thorne ^ in his analy- 
sis of the stability of neutron stars against radial collapse 
induced by an external tidal field. The motivation for 
Thome's paper was to refute the claims of Wilson, Math- 
ews, and Marronetti |Q regarding the "star-crushing" 
effect seen in their numerical simulations of binary neu- 
tron stars. We believe that our analysis has strengthened 
Thome's arguments. Ironically, Thorne and I p5[ now 
have reason, based on current quadrupole tidal coupling, 
to support the "star-crushing" effects observed by Wilson 
and Mathews fE6l in their revised simulations. 



This analysis has only been concerned with the tidal 
work involving a mass quadrupole moment interacting 
with an electric- type, quadrupolar tidal field in the slow 
motion approximation. It seems quite likely that this 
energy localization and gauge invariance of tidal work 
can be extended to higher order multipolar couplings. 
We will leave the demonstration of this for the future. 
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